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We study the nonequilibrium nanoscale energy transfer in a quantum junction in contact with fermionic and 
bosonic reservoirs, where the energy exchange is directly assisted by the electron-phonon interaction. The sys- 
tem describes the inelastic electron-phonon scattering near a nanoscale metal-dielectric interface and renders 
analytic investigation of the microscopic mechanism. We demonstrate the heat diode effect and negative dif- 
ferential thermal conductance in such nanoscale interfaces and discuss why these intriguing thermal properties 
are absent in bulk metal-dielectric interfaces. The results described here are relevant for understanding the in- 
terface heat transfer via electron-phonon coupling and could have potential applications in controlling energy in 
low-dimensional nanodevices. 
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In metal, electrons dominate heat conduction, whereas in 
dielectric phonons are main energy carriers [1, 2]. Hence, for 
metal-dielectric interfacial heat transport, energy transfer in- 
cludes three possible pathways: (1) energy exchange between 
phonons in metal and phonons in dielectric, which is widely 
studied with the acoustic and diffuse mismatch models [3-6] ; 

(2) nonequilibrium electron-phonon (e-ph) energy exchange 
within the metal, with subsequently phonon-phonon coupling 
across the interface [7-9]; (3) direct energy transfer from elec- 
trons in metal to phonons in dielectric through inelastic scat- 
tering induced by e-ph coupling across the interface [10-12]. 

When considering the pathway (2), one usually applies the 
phenomenological two-temperature model [13] that has been 
widely used to investigate the ultrafast pulsed laser heating 
of metals [14, 15], and assume the pathway (3) is negligible 
[8, 9]. However, when the imperfect interface or the large 
lattice-mismatch of metal and dielectric is present, both path- 
way (1) and (2) will be seriously suppressed and the pathway 

(3) becomes significant. Recent experiments have demon- 
strated that when studying the metal-dielectric interfacial heat 
transfer, the e-ph coupling across the interface becomes cru- 
cial [16, 17]. In particular, the experiment in Ref. [18] re- 
veals the strong direct coupling between electrons in metal 
and phonons in dielectric, which then dominates the overall 
interface heat transfer. A particular challenge that arises is 
then to resolve the microscopic mechanisms controlling the 
nanoscale interfacial energy transfer caused by the e-ph inter- 
action across the metal-dielectric interface. 

On the other hand, phononics recently emerges as a new 
discipline [19], where various functional thermal devices are 
designed to render the smart control of energy in micro-nano- 
scale. Among many intriguing properties, the heat diode ef- 
fect [20-25] and negative differential thermal conductance 
(NDTC) [23, 26-31] are the two most important ones for di- 
electric phononic devices acting as solid-state thermal rec- 
tifiers, switches, and transistors etc [19]. These dielectric 
phononic devices, when hybridized with metallic electronic 
circuits, then carry the potential of various beneficial applica- 
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FIG. 1. (Color online) (a) and (b) Sketch of inelastic scattering 
of electrons with losing (gaining) energy accompanied by phonon 
emission (absorption) due to the e-ph interaction across the metal- 
dielectric interface, (c) The two-level system setup describes the 
nanoscale metal-dielectric interface. 

tions, including the manipulation of heat dissipation and cool- 
ing in micro-nano-devices [32]. Therefore, understanding the 
interface heat transfer across nanoscale metal-dielectric hy- 
brid structures is a long-standing challenge not only of funda- 
mental interest, but also for practical applications in device 
design [6]. There are few previous efforts towards under- 
standing the interfacial heat transfer in bulk metal-dielectric 
systems through the direct e-ph interaction [10-12]. Never- 
theless, the heat diode effect and NDTC escape from the at- 
tention, and many important questions are left open: Can the 
bulk metal-dielectric with interfacial e-ph coupling possesses 
the heat diode and/or the NDTC properties? What if we scale 
the interface down to the micro-nano level? If so, what are 
their microscopic mechanisms in the nanoscale interfacial en- 
ergy transfer caused by the direct e-ph coupling? 

In this Letter, we shall answer these above mentioned ques- 
tions by studying a quantum junction in contact with elec- 
tronic and phononic reservoirs, where the energy exchange 
is directly assisted by the interfacial e-ph interaction, i.e. the 
pathway (3) [see Fig. 1]. Our simple model describes the es- 
sential inelastic e-ph scattering around metal-dielectric inter- 
faces in nanoscale and readily renders clear physical insights 
and analytic interpretations. Interesting phenomena like the 
heat diode effect and NDTC are identified in such nanoscale 
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metal-dielectric hybrid structures, which are overlooked in 
previous studies of macroscopic metal-dielectric interfacial 
heat transfer [10-12]. These effects are relevant for under- 
standing the interface heat transfer via e-ph interaction and 
would have potential applications in managing energy trans- 
fer in low-dimensional nanodevices. In the end, we explain 
the reason why the heat diode effect and NDTC are absent 
from the bulk metal-dielectric interfaces, and discuss under 
what conditions we can obtain these nontrivial properties. 

Generally, the e-ph interaction around the metal-dielectric 
interface is described by the Hamiltonian [1,2]: 



= ^M kq c\ +q c k (a q 

kq 



(1) 



where c k (c\) annihilates (creates) an electron of momentum 
k, a q (a q ) annihilates (creates) a phonon with momentum q, 
and Mk q is the e-ph coupling strength. This Hamiltonian de- 
picts two essential processes of inelastic e-ph scattering near 
the metal-dielectric interface that are responsible for the en- 
ergy transport: 1) an electron in the metal is down- scattered by 
emitting a phonon into the dielectric and the energy is trans- 
ferred from the metal to the dielectric [see Fig. 1(a)]; or 2) an 
electron in the metal is up-scattered by absorbing a phonon 
from the dielectric and the energy is transferred from the di- 
electric back to the metal [see Fig. 1(b)]. 

As the feature size of electronic and phononic devices con- 
tinues to scale down, as in the situation of quantum point 
contacts, the possible channels for the above scattering will 
be severely limited. In principle, when only considering a 
single-pair channel of the above two essential processes of the 
inelastic e-ph scattering, the nanoscale metal-dielectric inter- 
face can be described by a hybrid quantum junction with the 
Hamiltonian H = Hs + Hl + Hr + Vl + Vr, where a cen- 
tral system Hs interacts with two reservoirs H v (y = L,R) at 
different temperatures T v via the couplings V v [see Fig. 1(c)]. 
The central part is a two-level system: 



H s = stdldi 



E 2 d\d<. 



(2) 



where £\^ 2 denote two local energy levels with a gap ftwo = 
e 2 — S\. The left fermionic reservoir is described by the free 
electron Hamiltonian: 



(3) 



with the corresponding left coupling term describing the ex- 
change of electrons [33]: 

V L = J2 Xkl ( d i c fc+4 d i) +J2 Xk>2 ( d 2 c *' +4' d a) • (4) 



The right bosonic reservoir is modeled by independent 
phononic modes: 



(5) 



with the corresponding right coupling term describing the e- 
ph interaction that assists the relaxation and excitation of the 
central two-level system: 



V R 



(d\d 2 + d\di J (a q + at g J • 



(6) 



In the following, we consider the strong Coulomb- 
screening regime [34] so that the central system has only 
three eigenstates: empty state |0), state occupying the lower 
level |1), state occupying the upper level |2), with relations 
H s \0) = 0, H s \l) = ei|l>, H s \2) = e 2 \2). Employ- 
ing the quantum master equation in the interaction picture, 
the evolution of the total density matrix p nm = (n\p\m), 
(n, m G 0, 1, 2) is described by 



dp n 



dt 



imp(o)] r 



^fdr[V{tl[V{r)^{r)}} r 



(7) 

where V(t) is given in the interaction picture V(t) = 
e i(H s +H L +H R )t(y L + y R ) e -i(H s +H L +H R )t m Following the 

standard approach in open quantum systems [35], we arrive at 
the equation of motion for the three-state populations: 
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(8) 



as well as the heat flux from the left to the right: 

J = hu (p 2 k 2 ^i -2?ifci_> 2 ), (9) 

where p n = Tr^[p nn ] with Trg[-] tracing out the degrees of 
freedom of both baths and k n ^ m denotes the transition rate 
from state \n) to \m). 

The tunneling rates, describing the exchange 
of electrons between the left fermionic bath 
and the central system, are given by ko^ n = 

&Zk>ZnSZ, d ™-* nT/n nB[ct(T)c k (P)] = r L (e n )hn, 
kn^O = ^E^LIZo d ^ i£nT/H ^B[c k (T)c{(0)} = 

r L (e„)(l -f Ln ), where n= 1,2, f Ln = [ e ^-^)/(k B T L ) + 
l]" 1 and T L (e n ) = 2 fE k ^J(^n - e k ) = f \ 2 £n p L (e n ) 
with pl(s) the temperature-independent density of state 
(DOS) of the left electronic bath. 

The relaxation/excitation rates of the central system 
due to the inelastic e-ph coupling, accompanied by 
the phonon emission/absorption, are given by k 2 ^i = 
&E q MZS?„dTe™°TTr B [a q (T)al(Q)] = T R (u )(l + 
N R ),k^ 2 = ^E q M 2 q /r oo rfre--o-Tr B [at(r)a g (0)] = 
T r {uj )N r , where N R = [e^/^' - l]" 1 , r fi (w ) = 
X Z g M 2 S(fiw - ?u g ) = fMl o p R (hw ) with p R (hu) 
the temperature-independent DOS of the right phononic bath. 

By imposing the steady state condition po(t) = fti(t) = 
p 2 (t) = and the conservation law po(t) +pi(t) +p 2 (t) = 1 
to Eqs. (8) and (9), we finally obtain the heat flux: 

T h r r fL2 (1-/Ll) (1 + A^)~/l1 (1-/L2) N R 
J =ftOJ LLl-R ^ , 

(10) 
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where C = T R [(1 + f L2 ) (1 + N R ) + (1 + f L1 ) N R ] + 
(1 — /li/l2)- Note that if alternatively using the nonequi- 
librium Green's function method [37], we will arrive at the 
same results as above in the weak system-bath coupling limit. 

As we can see, the energy flux J is proportional to the dif- 
ference between two rate products: Jl2 (1 — + Nr) — 
hi (1 - /l 2 ) N r . The first product f L2 (1 - f L1 ) (1 + N R ) 
describes the relaxation rate of the central system from the 
occupied higher level to the empty lower one with releasing 
a phonon with energy Hujq into the right bosonic bath [see 
Fig. 1]. The second product /li (1 — /l2) Nr reversely de- 
scribes the excitation rate of the central system from the oc- 
cupied lower level to the empty higher one with absorbing a 
phonon with energy Hojq from the right bath. When Tl ^ Tr, 
the rate difference between these two processes gives rise to 
the nonequilibrium energy transfer. When Tl — Tr, these 
two processes are in balance and cancel with each other, so 
that the energy transport is absent. 

Set T L = T + AT/2, T R = T - AT/2, it can be 
verified from Eq. (10) that the thermal conductance k := 
dJ/dAT\^ T ^ , as a linear response quantity, is a single- 
valued function and always positive. However, as we will see 
in the nonlinear response regime, the heat diode and NDTC 
effects will emerge as intriguing transport properties of the hy- 
brid system. In fact, the heat flux can be generally expanded as 
J = ft AT + DAT 2 + O ( AT 3 ) , where the nonzero coefficient 
D as a result of the system asymmetry is responsible for the 
heat diode effect. Two intrinsic asymmetries are present in the 
system: i) two reservoirs aside the central system have differ- 
ent quantum statistics; ii) two system-bath couplings have dif- 
ferent natures, the left one is electron tunneling and the right 
one is e-ph coupling. Thus, considering these two fold asym- 
metries it is natural to have the asymmetric energy transfer 
when AT reverses its direction. 

When dJ/dAT = n + DAT + 0(AT 2 ) < 0, we say we 
have NDTC. The physical reason for the NDTC in our inter- 
face system is different from the previous mechanism in di- 
electric solid systems where phonon-phonon interaction (an- 
harmonicity, nonlinearality) is necessary to shift the phononic 
power spectrum [19] or in quantum junctions where the strong 
system-bath coupling is required [27]. Here, the system is un- 
der the weak system-bath coupling and the NDTC is a conse- 
quence of the competition between relaxation/excitation pro- 
cesses due to the e-ph coupling. In fact, from Eq. (10), we 
can verify that when fixing Tl, dJ/dAT = —OJ/OTr is al- 
ways positive, i.e., we can never have NDTC when we merely 
vary the temperature of the bosonic bath. For merely varying 
the fermionic bath temperature Tl, we find that dJ/dAT — 
OJ/OTl > under the condition e\ +ojq ^ \if ^ £i- There- 
fore, the NDTC is only possible when we at least vary the 
temperature of the left fermionic bath and when the two levels 
in the central are both either above or below the Fermi energy 
[ip of the left bath. 

Without loss of generality, in what follows we set the Fermi 
energy fip = as the energy reference point. Since spec- 
trum details of the DOS pL and pR of two reservoirs are not 
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FIG. 2. (Color online) Heat diode effect and NDTR in the two-level 
quantum junction, (a) Heat flux J as a function of the fermionic 
temperature Tl and the bosonic temperature Tr. Parameters are 
T L = F R = 2.5 meV/ft [36], si 25 meV, uj 50 meV. (b) 
The rectification ratio as a function of T for Tl = T ± AT/ 2, 
T R = T =F AT/2 with AT = 30 K. Inset depicts the rectifi- 
cation ratio as a function of the lower level energy at To = 200 
K and AT = 30 K. Other parameters are the same as above, 
(c) Heat flux vs AT along three corresponding lines in (a): (i) 
T L = To + AT, T R = T - AT; (ii) T L = T + AT, Tr = T ; 
(iii) T L = T ,Tr = T - AT, with T = 300 K. 



important for the present problem, the wide-band limit is ap- 
plied for both system-bath couplings so that T l and Tr are 
taken as constants. We plot a typical heat flux pattern as a 
function of both Tl and Tr in Fig. 2(a). It is seen that the 
heat flux pattern is asymmetric with respect to the separatrix 
Tl = Tr, which demonstrates clearly the heat diode effect. 
The rectification ratio \J {AT) / J {— AT) \ is defined to quan- 
tify the asymmetry of energy transfer under opposite temper- 
ature bias. As shown in Fig. 2(b), the heat diode effect be- 
comes significant when temperature is cooled down and the 
asymmetric heat transfer prefers the direction from electrons 
in metal to phonons in dielectric. In this way, the device acts 
like a good thermal conductor while in the reversed way, the 
device acts like a thermal insulator. The inset shows that the 
diode effect can be enhanced by tuning the level energy apart 
from the Fermi level /ip = 0. Moreover, it shows that the 
diode can be reversed near the Fermi level so that the asym- 
metric energy transfer prefers the direction from phonon in 
dielectric to electrons in metal, although this reverse effect is 
not distinct. Of course, enlarging the temperature difference 
can also enhance the diode effect due to its nonlinear response 
nature. 

For the behavior of another nonlinear response quantity, 
differential thermal conductance (DTC), Fig. 2(c) shows that 
for merely varying Tr [see the line (iii)], the system always 
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exhibits a positive DTC. Only with varying Tl [see the line 
(ii)] or both T L ,T R [see the line (i)], the NDTC emerges. 
This can be understood through Eq. (10): In the regime of 
Tr > Tl, when Tl approaches to the lower temperature 
regime, both /li,/l2 approach to if S\ + uoq > S\ > 
or /li, Il2 —> 1 if > S\ + c^o > e\. In this way, the relax- 
ation and excitation processes due to the e-ph interaction are 
effectively suppressed so that the energy flux \J\ deceases to 
zero, although the thermal bias Tr — Tl increases. 

In the regime of Tl > Tr, we can also have the NDTC 
effect if varying both temperatures. In fact, when Tl is large 
enough and Tr is low enough, we have /li ~ /l2 ~ 1/2 and 
N R * 0. In this way, J w l/^ r L iV(§r L + |r fl ). Now 
increasing Tr and the thermal bias Tl — Tr [e.g., see the line 
(iv) in Fig. 2(a)], /li, /l2 W1 U fix at the half value and TVr 
becomes some finite positive number. This modulation does 
not change the numerator of J but increases the denomina- 
tor, as: J « l^ r L r^/(|r L + |r fl + 3r^), so that 

I 



the energy flux deceases although the thermal bias increases. 
In fact, when Nr « 0, the relaxation and phonon emission 
is effective but the excitation and phonon absorption is sup- 
pressed. Thus, when Nr increases, the excitation becomes 
effective and absorbs the phonon from the right bosonic bath 
through the e-ph interaction which attenuates the energy trans- 
fer from metal to dielectric. Due to this physical reason, the 
NDTC appears in the regime of Tl > Tr as well if both tem- 
peratures are varied. 

Now let us examine why the heat diode and NDTC effects 
are overlooked in previous studies of the bulk metal-dielectric 
interface with direct e-ph coupling [10-12] and how can we 
observe them in the macroscopic interface. Starting form the 
Eq. (11), by using the Fermi golden rule [1, 10] or the stan- 
dard Bloch-Boltzmann-Peierls formulas [2, 38], one can get 
the energy transfer flux from the left electrons with Tl to the 
right phonons with Tr, as: 



poo poo 

J = 2ir dojF R (uj)hio / de PL {e) [f L (s + hw)(l - f L (e))(l + N r (Huj)) - f L (s)(l - f L (e + hw))N R (hLu)} , (11) 

JO J-oo 



where Pl(z) denotes the bulk electronic DOS and Fr(u) con- 
tains the bulk phononic DOS and the e-ph coupling matrix el- 
ements. This expression is familiar from earlier discussions of 
energy transfer between electron and phonon systems [1, 10- 
12, 38], and also contains essentially the same physics of the 
Eq. (10). To further proceed, one usually considers a good 
metal, i.e., the integral over de converges within a thermal 
energy fe^Tc around the Fermi level so that the bulk elec- 
tronic DOS can be taken as a constant pl(z) ~ Pl(pf)- 
Then, by applying the equality j defL^e + hu)(l — = 
hu/^W/ikBTL) _ = ^7v L (ftcj), one finally arrives 
[1, 11, 12, 38] at a Landauer-like formula, 

pOO 

J = 2k Pl (pf) / du)F R (u)(hu) 2 [N L (hcu) - N R (hu)} . 

J 

(12) 

Since the temperature dependence only manifests in the distri- 
bution difference Nl(Huj) — Nr(Huj), it is straightforward to 
verify that this bulk interfacial system can never have the heat 
diode and NDTC effects. This may also explain why the heat 
diode effect is absent in a bulk copper-copper oxide system 
[39]. 

From the above analysis, we see that the constant bulk 
DOS is a crucial assumption leading to the final Landauer- 
like formula. However, if pl(z) is not a smooth function but 
varies strongly on energy, which is exactly the case in low- 
dimensional nanoscale systems, we have to keep this elec- 
tronic DOS inside the energy integral. In this way, we will not 
get the Landauer-like Eq. (12) but instead retain the intrigu- 
ing properties of heat diode effect and NDTC. For example, 
assuming we have Pl(^) = Pl(pf)&(£ — Pf), then from Eq. 



(11), one arrives at 

poo 

J = 27tpl(pf) / dujT(uj) [N L (hw) - N R (hw)} ,(13) 

where T(oj) = F R (u)hLuk B T L In ^X^+i • As such ' 
the temperature dependence of T(uo) makes the heat diode 
and NDTC effects possible. In fact, the heat flux in Eq. (13) 
has qualitatively the same pattern of temperature dependences 
as in Fig. 2. Furthermore, note that the factor Fr(uj) contains 
the e-ph coupling matrix elements that are generally electronic 
energy dependent [1]. By considering non- smooth, strong 
energy dependent Fr(e,uj) in some designed interfacial sys- 
tems, we also expect to have the diode and NDTC effects. 

Finally, we note that our results can be extended to gen- 
eral fermionic-bosonic interfacial systems. For example, in 
the nanoscale metal-radiative field interfaces with electron- 
photon couplings, or in the nanoscale metal-magnetic insu- 
lator interfaces with electron-magnon couplings, we expect 
to observe the similar effects. These results described here 
are relevant for understanding the interface heat transfer and 
could have potential applications in controlling energy in low- 
dimensional nanode vices. 
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